Variable  Target  Model  Filtering  for  Use  with  Angle  Target  Tracking  Applications 


Philip  J.  Haney 
Raytheon  Syst^s  Ccwoipany 
Missile  Guidance  Lab 
Tewksbury,  MA.  01876 


Abstract 

A  variable  target  model  filter  is  investigated  for  use 
with  angle  target  tracking  applications.  In  the 
absence  of  a  target  maneuver,  a  "'quiescent**  target 
model  is  used  for  the  tracking  filter.  In  the  presence  of 
a  target  maneuver,  the  target  model  is  shifted  to  a 
"maneuvering**  model  The  basis  for  using  a 
"quiescent**  target  model  in  conjunction  with  a 
"maneuvering**  target  model  is  that  accurate  tracking 
performance  can  be  obtained  for  both  maneuvering 
and  non-maneuvering  situations  as  opposed  to  settling 
for  an  undesirable  compromise.  Simulation  results 
indicate  that  the  variable  target  model  filter  used  for 
this  work  is  capable  of  providing  "best  case**  miss- 
distance/tracking  performance  for  both  maneuvering 
and  non-maneuvering  target  conditions. 


Introduction 

A  distinctive  problem  with  many  tracking  filters  is  that 
their  target  models  are  quite  often  over-designed  for 
either  a  (near)  constant  velocity  target  or  a  (near) 
constant  acceleration  target.  As  a  result,  degradation 
in  performance  is  realized  when  a  tracking  filter 
containing  an  over-designed  target  model  is  applied  to 
a  class  of  targets  it  is  not  capable  of  handling.  For 
instance,  tracking  filths  strictly  designed  to  handle 
non-maneuvering  targets  will  suffer  in  the  presence  of 
a  target  maneuver  because  the  necessary  uncertainty 
(i.e.,  variance)  required  to  model  the  high^-order 
motion  of  the  target  is  not  available.  On  the  other 
hand,  tracking  filtos  strictly  designed  to  handle 
maneuvering  targets  will  suffer  in  the  presence  of  no 
target  maneuver  because  in  attempting  to  model  the 
(near)  zero-acceleration  component  of  the  target 
motion,  increased  estimation  errOTS  for  the  lower-order 
states  will  result. 


Consequently,  in  the  work  presented  hare,  the  target 
model  of  an  angle  target  tracking  filter  is  designed  to 
nominally  operate  in  a  “quiescent”  mode  during  the 
absence  of  a  target  maneuver.  However,  once  a  target 
maneuver  has  been  detected,  the  target  model  is  shifted 
to  a  “maneuvering”  model.  The  premise  for  switching 
between  a  “quiescent”  target  model  and  a 
“maneuvering”  target  model  is  that  accurate  tracking 
performance  can  be  obtained  for  both  situations  as 
opposed  to  settling  for  an  undesirable  compromise. 
Simulation  results  illustrate  that  the  angle-frame 
variable  target  model  Iracking  filto*  used  for  this  work 
is  capable  of  providing  “best-case”  miss- 
distance/tracking  performance  (i.e.,  “best-case”  with 
respect  to  a  target  model  optimized  for  a  specific 
maneuver  condition)  for  both  maneuvering  and  non¬ 
maneuvering  target  conditions. 

The  breakdown  of  this  work  is  outlined  as  follows: 
The  standard  target  equations  of  motion  (continuous¬ 
time  and  discrete-time)  which  will  be  refeenced 
throughout  this  paper  are  provided  in  the  opening 
section.  Next,  the  discrete-time  Kalman  filter 
estimator  equations  used  for  modeling  the  standard 
target  equations  of  motion  are  given.  The  “quiescent” 
target  model,  the  “maneuvering”  target  model  and  the 
target  detection  scheme  utilized  for  this  work  are  then 
discussed  in  subsequent  sections,  and  finally, 
simulation  results  conclude  the  papa*. 


Target  Motion  State  Models 

The  standard,  continuous-time  target  equation  of 
motion  is  described  as  follows*: 

jt(0  =  F{t)x(t)  +  Gw{t)  (1) 


For  sin^licity,  a  single  angle-frame  tracking  filter  is  referred  to 
throughout  this  discussion.  However,  individual  azimuth  and  elevation 

_  angle-frame  tracking  filters  would  most  likely  be  desired  for  actual 
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Extending  equation  1  to  disaete-time  yields  the 
following: 

x{k  + 1)  =  +  u(Jc)  (2) 

whCTC  O  is  the  target  state  transition  matrix  and  u  is 
an  inhomogeneous  driving  input.  It  should  be  noted 
that  even  though  the  input  u  is  not  a  sampled  version 
of  the  continuous-time  white  noise  input  w  of 
equation  1,  it  can  be  shown  that  u  is  white  in  the 
discrete  time  sense.  Consequently,  equation  2  is 
directly  suitable  for  Kalman  filter  applications.^ 
Furthermore,  since 

r+Af 

t 

it  can  be  seen  that  for  the  disaete-time  target  state 
model  provided  in  equaticm  2 

(4) 

(i5:+l)A/ 

tt(fc)=  (5) 

kAt 


Kalman  Filter  Equations 

Fa  this  work,  it  is  assumed  that  an  active  tracking 
sensor  provides  target  position  based  upon  the 
following  angle  measurement: 

y{k)  =  Hx(k)  +  v(k)  (6) 

where  H  is  the  observation  matrix  and  v  is  additive 
white  noise  with  variance  R  . 

Equations  2  and  6  provide  the  necessary  plant  and 
observation  models  for  the  implementation  of  a 
discrete-time  Kalman  filter.  The  Kalman  filter 
estimator  equations  used  for  this  work  utilize  the 
following  notation:'* 

1.  Correction  (based  upon  new  measurement): 

(7) 

P^i+)  =  [I-K^H^]P,(-)  (9) 


2.  Prnpagatirat  (to  next  measurement) 
h+i  (~)  =  (+) 

Pk^i(-)=^kPkW%^  +  Qk  (11) 

where  K  is  the  q)timal  Kalman  gain  matrix,  P  is  the 
state  estimate  error  covariance  matrix,  y,  H  and  R 
are  taken  from  equation  6,  4>  is  taken  from  equation 
4,  and  Q  is  the  process  noise  covariance  matrix  based 
upon  equation  5  and  defined  as  follows:^ 

e, =£{«,«/)  (12) 

Equations  7-12  provide  the  angle-frame  target  tracking 
state  estimator  equations  used  in  this  work  for  both  the 
‘‘quiescent”  and  “maneuvering”  target  models. 

Quiescent  Target  Model 

In  the  absence  of  a  target  maneuver,  the  target  model  is 
placed  in  a  “quiescent”  mode  of  operation.  The 
“quiescent”  target  acceleration  is  modeled  using  a  first- 
order  markov  process  and  process  noise  to  allow  for 
slight  changes  in  the  target’s  motion.  The  “quiescent” 
target  model  is  defined  by  the  following  continuous¬ 
time  parameters  based  upon  equation  1:^’^ 
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=  quiescent  white  noise  driving  function  with 
power  spectrum - — 

P*Tq 

cTjy  =  quiescent  target  accelaation  standard  deviation 

=  quiescent  target  acceleration  time  constant 

=  inverse  of  quiescent  target  acceleration  time 
constant 

r  =  relative  range  measurement/estimate 

Considering  equations  4,  5,  and  12,  the  “quiescent” 
target  model  parameters  required  fa  the 
implementation  of  a  disaete-time  Kalman  filta  are 

given  as  follows  (assuming  Af  «  ): 
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r  =  relative  range  measurement/estiinate 

Again, .  considering  equations  4,  5,  and  12,  the 
“maneuvering”  target  model  parametCTS  required  for 
the  implementation  of  a  discrete-time  Kalman  filt^  are 
given  as  follows  (assuming  Ai  «  ): 
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Maneuvering  Target  Model 


When  a  target  maneuver  has  been  detected,  the  target 
model  is  shifted  to  a  “maneuvering”  mode  of 
operati(Mi.  The  “maneuvering”  target  acceleration  is 
modeled  using  a  first-OTder  markov  process  and 
process  noise  to  allow  for  severe  changes  in  the 
target’s  motion.  The  “maneuvering”  target  model  is 
defined  by  the  following  continuous-time  parameters 
based  upon  equation  1:^’^ 
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=  maneuver  white  noise  driving  fimction  with 
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power  spectrum 


=  maneuver  target  acceleration  standard 
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=  maneuver  target  acceleration  time  constant 

=  inverse  of  maneuver  target  acceleration  time 
constant 
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Target  Maneuver  Detection 

The  target  maneuver  detection  scheme  used  for  this 
work  is  based  on  a  fading  memory  average  of  the 
“quiescent”  target  model’s  normalized  innovations  as 
follows: 

[y{k)-xAk)f 

where:  y  =  line-of-sight  measurement 
jcj  =  line-of-sight  state  estimate 
Pn  =  line-of-sight  state  estimate  error 
covariance 

R  =  line-of-sight  measurement  noise 
variance  estimate 

0  <71,72  <1 
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Basically,  assuming  a  “quiescent”  mode  of  operation, 
on  the  occasicHi  that  equaticm  13  exceeds  a  given 
threshold  it  is  accepted  that  a  target  maneuver  is 
present,  and  the  “maneuv^ing”  target  model  is 
transitioned  to  from  the  “quiescent”  target  model.  The 
“quiescent”  target  model  is  th^  returned  to  if  the 
estimated  target  accelerations  over  a  certain  “window” 
of  time  become  statistically  insignificant  with  respect 
to  their  standard  deviations.  The  test  statistic  used  to 
determine  the  significance  of  the  target  acceleration 
estimates  while  in  a  “maneuvering”  mode  of  operation 
is  given  as  follows: 

k 

p(k)=  (14) 

j=k-l+l 


Simulation  Results 

In  order  to  evaluate  the  performance  of  the  angle-frame 
variable  target  model  filter,  various  parameters  (i.e., 

required  for  the  “quiescent”  and 

“maneuvering”  modes  of  operation  need  to  be 
determined. 


where:  S(j) 

P 

33 
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^  NjU)^ 

^33  0) 

=  target  acceleration  state 
estimate 

=  target  acceleratitHi  state 

estimate  error  covariance 
=  window  length 


Consequently,  for  a  specific  window  length,  I ,  when 
the  sum  of  equation  14  exceeds  unity,  the  target 
acceleration  is  deemed  significant  with  respect  to  its 
standard  deviation  and  vice-versa.  Furthermore,  in 
order  to  allow  for  a  smoother  transition  between  the 
“quiescent”  and  “maneuvering”  target  models,  once 
the  “maneuvering”  target  model  has  been  initiated, 
each  time  the  resultant  summation  of  equation  14 
exceeds  unity  a  maneuver  probability  value  is 
incremented  by  0.1.  Similarly,  each  time  the  resultant 
summation  of  equation  14  is  less  than  unity  the  same 
maneuver  probability  value  is  decremented  by  0.1, 
Basically,  this  maneuver  probability  value  represents 
the  severity  of  the  target  maneuver  and  can  be  used  to 
weight  the  “maneuvering”  target  model’s  process  noise 
in  order  to  provide  a  more  seamless  transition  between 
target  models.  In  addition,  once  the  maneuver 
probability  reaches  zero  the  “quiescent”  target  model  is 
reinstated. 


The  “maneuvering”  target  model’s  process  noise 
weighted  by  the  aforementioned  maneuver  probability 
is  given  as  follows: 


Consequently,  for  the  “quiescent”  target  model,  if  the 
normal  target  accel^ation  (Nj)  is  assumed  to  be  a 
random  variable  uniformly  distributed  over  the  interval 
(-Nr  ),  then  its  corresponding  variance  can 

be  computed  as  follows:^ 


IN^ 
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^dx~- 


N-r 


(15) 


-Nrr 


During  a  “quiescent”  mode  of  operation,  any  change  in 
target  motion  is  assumed  to  be  solely  attributed  to 
atmospheric  turbulence/wind  effects.  Consequently, 
the  maximum  expected  “quiescent”  lateral  target 
acceleration  used  for  this  work  is  assumed  to  be  1  gee. 
As  ^  .wouit,  the  target  acceleration  variance  required 
for  computing  the  “quiescent”  target  model’s  process 
noise  is  given  as  follows: 


(1.0-9.807/ 

3 


=  32.059 


(m^/sec"^) 


For  a  “maneuvering”  mode  of  operation,  any  change  in 
target  motion  is  believed  to  be  the  result  of  an 
intentional  target  maneuver.  For  this  work,  the 
maximum  expected  “maneuvOTng”  lateral  target 
acceleration  is  assumed  to  be ’6  gees.  Consequently, 
the  target  acceleration  variance  required  for  computing 
the  “maneuvering”  target  model’s  process  noise  is 
given  as  follows: 


2 


(6.0-9.807/ 
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1154.1270 


(m^/sec"^) 
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The  time  constants  used  for  this  work  refer  to 

the  duration  of  the  expected  “target  maneuver”.  For  a 
“quiescent”  mode  of  operation,  typical  time  constants 
corresponding  to  atmospheric  turbuience/wind  effects 
have  been  shown  to  be  on  the  order  of  1-2  seconds.^ 
Consequently,  for  this  work  a  time  constant  of  2 
seconds  (i.e.,  =  0.5)  is  used  to  model  the 

atmospheric  turbulence/wind  effects  for  the 
“quiescent”  target  model.  Similarly,  for  a 

“maneuvering”  mode  of  operation,  typical  target 
maneuver  time  constants  for  the  application  under 
consideration  can  range  anywhere  from  1-10  seconds. 
As  a  result,  for  this  work  a  time  constant  of  4  seconds 
(i.e.,  =  0.25)  is  used  to  model  the  intentional 

target  maneuver  eftects  for  the  “maneuvering”  target 
model. 

The  simulation  results  presented  in  the  following 
figures  compare  Mcmte  Carlo  performance  data 
obtained  for  a  constant  “quiescent”  target  model,  a 
constant  “maneuvering”  target  model,  and  a  variable 
target  model  corresponding  to  two  different  target 
maneuver  scenarios.  The  two  target  maneuver 
scenarios  considered  for  this  work  are  described  as 
follows: 

1.  No  target  maneuver. 

2.  Evasive  elevation  maneuver  (6gee  pullup)  induced 
at  a  time-to-go  of  3  seconds. 

In  addition  to  the  above,  the  following  target  initial 
conditions  and  parameters  further  describe  the 
simulation  scenarios  investigated  for  this  work: 

1 .  Target  initial  downrange  of  25km. 

2.  Target  initial  altitude  of  lOOm. 

3.  Target  velocity  of  150  m/s. 

4.  Target  flying  straight  in  along  primary  target  line 
(PTL). 

Figure  *1  illustrates  the  miss-distance  performance  of 
the  various  target  models  corresponding  to  target 
maneuver  scenario  ^1.  From  this  figure  it  can  be  seen 
that  all  three  target  models  essentially  provide  similar 
miss-distance  performance. 


Figure  ^1  Normalized  Miss-Distance  Performance 
For  Target  Maneuver  Scenario  ^1 


Table  I  provides  the  tracking  performance  of  the 
various  target  models  corresponding  to  target 
maneuver  scenario  ^1.  The  performance  index,  7,  used 
for  this  work  is  described  as  follows: 

N  I 

^\xik)-xXk) 

J  ^  ^ ^ -  (16) 

'■  N 

where:  X;  =  i"'  true  state 

X-  =  i"'  state  estimate 
N  =  number  of  data  samples 
i  =  state  number  (i  =  1,2,3) 

In  addition  to  the  individual  tracking  performance 
indices  of  each  state  and  each  filter.  Table  I  also 
provides  an  overall  tracking  performance  index  which 
is  calculated  as  the  sum  of  the  individual  azimuth  and 
elevation  tracking  performance  indices  corresponding 
to  each  of  the  three  states  within  the  tracking  filter. 
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Table  I  Tracking  Performance  For 
Target  Maneuver  Scenario 


Constant- 

Quiescent 

Model 

Constant- 

ManeuvOT 

Model 

Variable 

Model 

LOS 

Azimuth 

Elevation 

3.4218e-5 

L3610e-4 

4.7793e-5 

1.3869e-4 

3.7204e-5 

1.3822^ 

LOS  Rate 

Azimuth 

Elevation 

1.4895e-4 

2.6890&4 

2.5301e-4 
2.7739^  : 

L3937e-4 

2.4813e-4 

Acceleration 

Aaimuth 

Elevation 

0.0910 

0.0652 

0.7410 

0.4016 

0.0976 

0.0961 

Overall 

Perfbnnance 
(Sum  of  J/s) 

0.1568 

1.1433 

0.1943 

From  Table  I  it  can  be  seen  that  the  tracking 
performance  of  the  “constant-maneuver”  target  model 
suffers  in  the  absence  of  a  target  maneuver  while  the 
“constant-quiescent”  target  model  and  the  “variable” 
target  model  track  the  non-maneuvering  target  quite 
well. 

Finally,  Rgure  *2  illustrates  the  line-of-sight  rate 
tracking  performance  of  the  various  target  models 
corresponding  to  target  maneuver  scenario  ^1.  From 
this  figure  it  can  be  seen  that  the  line-of-sight  rate 
estimates  provided  by  the  “constant-maneuver”  target 
model  are  fairly  “noisy”  while  those  provided  by  the 
“constant-quiescent”  target  model  and  the  “variable” 
target  model  are  quite  clean.  The  reason  for  the 
“noise”  observed  on  the  “constant-maneuver”  target 
model’s  line-of-sight  rate  estimates  is  that  while 
attempting  tc  "t  ’  ’  the  (near)  zero  acceleration 
component  of  the  target  motion,  the  high-level  of 
target  maneuver  uncCTtainty  inherent  within  the 
“constant-maneuver”  target  model  induces  increased 
estimation  errors  for  the  lower-order  states. 

Figure  *3  displays  the  miss-distance  performance  of  the 
various  target  models  corresponding  to  target 
maneuver  scenario  ^2.  From  this  figure  it  can  be  seen 
that  the  “constant-quiescent”  target  model  significantly 
suffers  in  the  presence  of  a  target  maneuver  while  the 
“constant-maneuver”  target  model  and  the  “variable” 
target  model  perform  quite  well. 


Constant  Quiescail  Target  Modd 


VallW^ 


Figure  ^2  LOS  Rate  Tracking  Performance 
For  Target  Maneuver  Scenario  ^1 
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Figure  ^3  Normalized  Miss-Distance  PerfMmance 
For  Target  Maneuver  Scenario  *2 
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Table  11  provides  the  tracking  performance  of  the 
various  target  models  corresponding  to  target 
maneuver  scenario 


Table  n  Tracking  Performance  For 
Target  Maneuver  Scenario  *2 


Constant- 

C^iiescent 

Model 

Constant- 

Maneuv^ 

Model 

Variable 

Model 

LOS 

Azimuth 

Elevation 

5.5317e-5 

1.4245e4 

4,9445e-5 

1.3285e4 

4.5561e-5 

1.3008e-4 

LOS  Rate 

Azimuth 

Elevation 

2.5661e-f 

9.7831e-4 

2.3578e4 

3.8342e4 

2.0963e4 

3.7928e4 

Accel^ation 

Azimuth 

Elevation 

0.1242 

1.1370 

0.6918 

0.7810 

0.1410 

0.6826 

Overall 

Performance 
(Sum  of  Ji’s) 

1.2626 

1.4736 

0.8244 

From  Table  II  it  appears  that  the  “constant-quiescent” 
target  model  tracks  the  maneuvering  target  better  than 
the  “constant-maneuver”  target  model.  However,  a 
closer  look  reveals  an  explanation  for  the  seemingl)^ 
misleading  results.  Basically,  prior  to  the  end-game 
maneuver,  the  “constant-quiescent”  target  model 
provides  the  better  tracking  performance  in  both 
channels  (i.e.,  azimuth  and  elevation)  while  during  the 
target  maneuver,  the  “constant-maneuver”  target 
model  provides  the  better  tracking  performance  i  the 
maneuver  channel  only.  Consequently,  when 
considering  both  channels  throughout  the  duration  of 
the  flight,  the  tracking  performance  of  the  “constarit- 
quiescent”  and  “constant-maneuver”  target  models 
becomes  “compromised”  which  leads  to  similar 
“overall”  tracking  performance  indices.  On  the  other 
hand,  it  can  be  seen  from  Table  II  that  the  “variable” 
target  model  is  capable  of  providing  “best-case” 
tracking  performance  for  both  non-maneuvering  and 
maneuvering  target  conditions.  For  instance,  prior  to 
the  end-game  maneuver,  the  “variable”  target  model  is 
capable  of  providing  tracking  performance  comparable 
to  the  “constant-quiescent”  target  model  while  during 
the  target  maneuver,  the  “variable”  target  model  is 
capable  of  providing  tracking  performance  comparable 
to  the  “constant-maneuver”  target  model  in  the 
maneuver  channel  and  the  “constant-quiescent”  target 
model  in  the  non-maneuver  channel. 


Finally,  Figure  ^4  displays  the  line-of-sight  rate 
tracking  p^formance  of  the  various  target  models 
correspcHiding  to  target  maneuver  scenario  ^2.  From 
this  figure  it  is  observed  that,  again,  the  “constant- 
quiescent”  target  model  provides  the  better  tracking 
performance  prior  to  the  end-game  maneuver  while  the 
“constant-maneuver”  target  model  provides  the  better 
tracking  performance  during  the  maneuver  (in  the 
maneuver  channel  only).  Furthermore,  it  can  also  be 
seen  that  the  “variable”  target  model  is  capable  of 
providing  “best-case”  tracking  performance  (in  both 
channels)  prior  to  and  during  the  end-game  maneuver. 


()](ist2iit(tertTar9£tM^ 


Figure  *4  LOS  Rate  Tracking  Performance 
For  Target  Maneuver  Scenario  *2 
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Conclusion 


A  variable  target  model  Gita'  has  be^  investigated  for 
use  with  angle  target  tracking  applications.  The 
filt»'’s  target  model  is  shifted  between  a  “quiescent” 
mode  of  operation  and  a  “maneuvering”  mode  of 
operatimi  based  upon  the  absraice  or  presence  of  a 
target  maneuver.  The  basis  for  using  a  “quiescent” 
target  model  in  exjunction  with  a  “maneuvering” 
target  model  is  that  accurate  tracking  perfiamance  can 
be  obtained  for  both  maneuvering  and  non¬ 
maneuvering  situations  as  opposed  to  settling  for  an 
undesirable  compreanise.  Simulation  results  illustrate 
that  the  variable  target  model  filto"  is  capable  of 
providing  “best-case”  miss-distance/tracking 
perfcamance  (i.e.,  “best-case”  with  respect  to  a  target 
model  q)timized  for  a  specific  maneuvCT  cxdition)  for 
both  maneuvering  and  nx-maneuvering  target 
xnditions. 
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